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Abstract. We consider Potts model hypersurfaces defined by the multivariate 
Tutte polynomial of graphs (Potts model partition function). We focus on the 
behavior of the number of points over finite fields for these hypersurfaces, in 
comparison with the graph hypersurfaces of perturbative quantum field theory 
defined by the Kirchhoff graph polynomial. We give a very simple example of 
the failure of the "fibration condition" in the dependence of the Grothendicck 
class on the number of spin states and of the polynomial countability condition 
for these Potts model hypersurfaces. We then show that a period computation, 
formally similar to the parametric Feynman integrals of quantum field theory, 
arises by considering certain thermodynamic averages. One can show that 
these evaluate to combinations of multiple zeta values for Potts models on 
polygon polymer chains, while silicate tetrahedral chains provide a candidate 
for a possible occurrence of non-mixed Tate periods. 



1. Introduction 

A lot of attention was devoted in recent years to the intriguing occurrence of 
periods of algebraic varieties and motives in the context of perturbative quantum 
field theory and residues of Feynman integrals. The question initially arose from 
numerical observations of Broadhurst and Kreimer |12j on the occurrence of mul- 
tiple zeta values in the computation of Feynman diagrams of massless scalar field 
theories. Multiple zeta values were conjectured to be the periods of mixed Tate 
motives, a result recently proved by Francis Brown |13j . Moreover, the parametric 
formulation of Feynman integral (see [9], [53]) exhibits the residues of Feynman 
diagrams explicitly as periods of an algebraic variety (the graph hypersurface com- 
plement), at least modulo divergences and renormalization. It was then natural to 
expect that the occurrence of multiple zeta values would be an indication of the 
fact that the graph hypersurfaces are mixed Tate motives. This was conjectured by 
Kontsevich, in a substantially equivalent statement predicting that these hypersur- 
faces would be polynomially countable. This means that their counting of points 
over finite fields would depend polynomially on the order of the field. At first this 
conjecture was confirmed (see [35], [3S]) for all graphs with up to twelve edges, but 
it was later disproved by a very elegant and general result of Belkale and Brosnan 
[7], which shows that, contrary to expectations, the classes of the graph hypersur- 
faces span a suitable localization of the Grothendieck ring of varieties, hence they 
can be arbitrarily far from the mixed Tate case as motives. For a discussion of the 
difference between the behavior of the classes in the localized or in the ordinary 
Grothendieck ring see also [2] . The original cases observed numerically in [T2] were 
later confirmed to be periods of mixed Tate motives in [14], [15], [33] and other 
cases in which the parametric Feynman integral turns out to be a period of a mixed 
Tate motive were analyzed in [3], [IS], [H], [IE], [31], [32], and others. 
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While the result of Bclkalc and Brosnan |i7j that disproved the Kontsevich con- 
jecture is very general and elegant, it is not immediately constructive, in the sense 
that it does not exhibit an explicit graph that fails to satisfy the polynomial count- 
ability property. The first such explicit example was identified by Doryn in |17j . 
and other examples were then constructed in |16j . There are interesting conceptual 
differences between these examples, which point to two different forms in which the 
polynomial countability property can fail. One, as in the case of |17j . where the 
dependence in the order of the field fails to satisfy the same polynomial function at 
all primes, and another more drastic one, where an actual non-polynomial function 
appears. 

In the recent paper pQ, it was shown that some of the techniques employed in 
the study of the classes of the graph hypersurfaces in the Grothendieck ring of 
varieties can be extended to another type of physical model, this time not in quan- 
tum field theory but in statistical mechanics. In fact, the partition function for 
the Potts model for a spin system on a graph with q possible states is given by a 
combinatorial polynomial associated to the graph, the multivariate Tutte polyno- 
mial, see [31]. This is closely related to the Kirchhoff polynomial that defines the 
graph hypersurfaces in quantum field theory. In the context of Potts models, one 
is less directly interested in the periods of the hypersurface complement, and more 
in describing how the hypersurface behaves when over families of graphs that ap- 
proximate some infinite graph. In fact, zeros of the partition function correspond 
to physical phase transitions, if they happen for non-negative real values of the 
parameters, while they can be thought of as some kind of non-physical "virtual" 
phase transitions otherwise. While (at least in the ferromagnetic case) the Potts 
model on a finite graph does not have phase transitions, which can only appear in 
the limit of an infinite graph, in the antiferromagnetic case there are possible phase 
transitions also over a finite graph [31] . The sets of real zeros of the multivariate 
Tutte polynomials have been extensively studied in the Potts model literature, most 
notably by Jackson and Sokal [53] ■ In pQ it is shown that one can estimate how 
the complexity of this set of zeros grows in certain families of graphs that approach 
infinite graphs, in terms of classes in the Grothendieck ring of the varieties defined 
by the multivariate Tutte polynomials, the Potts model hypersurfaces. These are 
analyzed using a deletion-contraction type formula similar to the one obtained in 
[3] for the graph hypersurfaces of quantum field theory. 

In the present paper, continuing to draw on the analogy between the graph 
hypersurfaces and the Potts model hypersurfaces, Zc, q , for a graph G and for a 
number q of spin states, we compare the behavior of the counting of points over finite 
fields for these varieties. We first consider the "fibration condition" discussed in [T], 
that is, the question of whether the Grothendieck class of Za,q depends on q outside 
of the special values q = and q = l. This question is motivated by an observation 
made in pQ, where it was shown that, for certain classes of graphs (such as polymer 
chains), the classes of the hypersurfaces Z^ q behave as if the hypersurfaces were 
a locally trivial fibration away from those two special values. We show here that 
the fibration condition already fails for the tetrahedron graph KA, and we argue 
that one can use a Monte Carlo method, similar to the one developed in [57], 
[23] . to test the failure of the fibration condition for Potts model hypersurfaces of 
more complicated graphs. We then look at the polynomial countability question for 
the Potts model hypersurfaces. Naturally, one expects that failures of polynomial 
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countability will appear for much smaller and combinatorially simpler graphs than 
in the cases of the graph hypersurfaces of [T5] and [T7]. Indeed, this is the case. A 
trivial kind of failure of polynomial countability comes from the primes that divide 
the number of spin states q, but even after reformulating the question by taking into 
account this trivial failure, we see that, again, the tetrahedron graph KA already 
provides an example where a more interesting failure of polynomial countability is 
possible. 

We discuss the role of periods in this statistical mechanical setting, where we in- 
terpret them in terms of certain thermodynamic averages. Unlike the more difficult 
case of quantum field theories, where infrared divergences arise from the intersection 
between the locus of integration in the parameteric Feynman integral and the graph 
hypersurfaces, in the (ferromagnetic) Potts model case the hypersurface does not 
intersect the simplex over which the integration is performed, hence we can directly 
interpret the resulting integral as a period, without the serious additional compli- 
cation of blowups performing regularization and renormalization, as one finds ( |10j , 
|11| ) in the quantum field theory case. 

The main purpose of this paper is pedagogical. It is aimed at presenting this 
approach to Potts models through the algebro-geometric and motivic properties of 
the hypersurfaces defined by the multivariate Tutte polynomials, in terms of very 
concrete and simple examples, which, we hope, can serve as an illustration of the 
general methodology. 

2. Graph hypersurfaces and Potts model hypersurfaces 

To fix notation, in the following we denote by G = G(V, E) a finite graph, with 
vertex set V and edge set E. We also write V = V(G) and E = E(G). 

2.1. Quantum field theory and graph hypersurfaces. We first recall briefly 
the definition of the graph hypersurfaces of perturbative quantum field theory. 

Definition 2.1. Let G = G(V, E) be a graph with n — #V . Assign an edge weight 
(a variable) t e to each edge e G E. The Kirchhoff graph polynomial of G is of the 
form 



where the sum is over all the spanning trees T (or maximal spanning forests in 
the multi- connected case). The graph hypersurface Xq is the hypersurface in A n 
defined by the vanishing of the graph polynomial 



Since the graph polynomial is homogeneous of degree 6i(G), the affine graph 
hypersurface is the affine cone over a projective graph hypersurface in P" . 

The role of the graph hypersurfaces in perturbative quantum field theory is ex- 
plained by the following result (see [§], [21], The Feynman integral associated 
to the graph G is given by 



(2.1) 




TcG e^E(T) 



(2.2) 



X G = {t = (ii,..., t n ) G A" | * G (t) = 0}. 



(2.3) 
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where D is the spacetime dimension, Pq is the second graph polynomial (whose 
explicit form we do not recall here), with p the external momenta, a n the unit 
simplex, and oj n the volume form. Here the Feynman integral refers to the residue 
(after removing a divergent Gamma factor), written in momentum space and in 
the Feynman parametric form, for the case of a massless scalar field theory (see the 
references listed above for more details). Modulo important issues of divergence 
and renormalization, the integral (2.3) is then a period of the graph hypersurface 
complement. 



2.2. Statistical mechanics and Potts model hypersurfaces. The analogous 
object we will be discussing in the statistical mechanics setting is the Potts model 
partition function. 

To avoid notational ambiguity, we will write p for an integer prime and q = p n 
for a prime power, so that we will be using the notation F q is the finite field with 
q elements. We reserve the notation q for the variable in the Potts model partition 
function that denotes the number of spin states. 

Definition 2.2. Let G = G(V, E) be a graph with n = #V , and with edge weights 
t = (t e )eev- The multivariate Tutte polynomial of G is 

(2.4) T G (q,t)= 5> fe(A) n*- 

ACE eSA 

where the sum is over the subgraphs of G that span all vertices ( that is, subgraphs 
G' with V(G') = V{G) and E(G') = A C E), and where k{A) is the number of 
connected components in each subgraph, and q is an indeterminate. 

To see the relation to Potts models in statistical mechanics, we consider a set 21 
of cardinality q (here assumed to be a non-negative integer) of possible spin states, 
and a statistical system of spins assigned to the vertices of a graph G. A state of 
the spin system if an assignment of a spin state a v to each vertex. The energy of a 
state is the sum over all edges of a quantity 



H(e) = 



o~ v ^ o~ w 



where 9(e) = {v, w}, and where J e is a fixed value, with J e > in the ferromagnetic 
case and — oo < J e < in the antiferromagnetic case. The weight edges are then 
related to the energies J e by the relation 

t e = eP J > 1, 

where f3 is (up to the Boltzmann constant) an inverse temperature. The partition 
function of the Potts model statistical system onG = G(V, E) is then given by the 
expression 

(2-5) Z G {q,t)= ]T H(l + t e 6^ <<T J, 

cr:V-J-2l eEE 

where the sum is over all states (that is, all maps of vertices to spin states), 5 is 
the Kronecker delta, and 9(e) = {v,w}. The relation of the multivariate Tutte 
polynomial ( |2.4[ ) to statistical mechanics comes from a famous result of Fortuin- 
Kasteleyn [20 , which shows that the partition function of the Potts model is a 
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restriction of the multivariate Tutte polynomial 

(2-6) Z G (q,t) =T G (q,t)\ qen>t=e 0j_ leR . 

The Potts model hypersurface is the locus of zeros defined by the multivariate 
Tutte polynomial 

(2.7) Z G = {(q,t)G A n+1 \Z G (q,t) = 0}. 

Notice that, unlike the Kirchhoff graph polynomial, the multivariate Tutte poly- 
nomial is not homogeneous, so this is an affine hypersurface in A n+ , where n = 
#£'(G), but not a projective hypersurface. One also considers the Potts model 
hypersurface for fixed q, 

(2.8) Z G>q = {t€A n \Z G (q,t)=Q}. 

The relation between the Potts model hypersurfaces and the graph hypersurfaces 
of quantum field theory lies in the fact (see pQ, [M]) that the Kirchhoff graph 
polynomial ^ G , or rather its equivalent form 

(2.9) * G (t) =j2 n 

TcG e£E(T) 

which is related to ^ G by dividing by IlegB ^ e an< ^ a PPly m g the Cremona trans- 
formation t e i — ^ l/t e to all the edge variables, is obtained by normalizing the 
T G (q,t) = q~ k ( G )T G (q,t), evaluating at q — and taking the homogeneous piece 
of lowest degree. This relation has a more geometric interpretation in terms of 
tangent cones, as in Lemma 2.7 of pQ. This relation is the motivation for extending 
techniques developed in the quantum field theory context for studying the geometry 
of the graph hypersurfaces to the setting of Potts models. 

2.3. Counting points and classes in the Grothendieck ring. Let Vk de- 
note the category of (quasi-projective) algebraic varieties over a field K. The 
Grothendieck ring of varieties Kq(Vk) is the free abelian group generated by the 
isomorphic classes [X] of varieties, modulo the relations 

[X] = [Y] + [X x Y] 

for closed embeddings Y C X of subvarieties. The product operation that gives the 
ring structure is given by [X] ■ [Y] = [X x Y] . 

For an affine hypersurface X C A™, we will use in the following the notation 
{X} for the class in -Ko(Vk) of the hypersurface complement 

(2.10) {X} := [A™ \X} = [A"] - [X] = L" - [X], 
where L = [A 1 ] is the Lefschetz motive, the class of the affine line. 

Similarly, one can define a Grothendieck ring Kq{Vz) for arithmetic varieties 
defined over Z. While the relations are formally the same, the classes [X] in this 
case denote isomorphism classes as varieties over Z. We will sometime write only 
Kq(V), without specifying the field or ring of definition, when it should be clear 
from the context. 

An additive invariant is a map x '■ Vk - > R-, with values in a commutative ring R, 
that satisfies x(X) = xOO if X and Y are isomorphic, and x(X) = x(^)+x(^ x ^) 5 
for closed embeddings Y C X. Moreover, one requires that x{X xF) = x(Jf)x(y). 
An additive invariant determines and is determined by a ring homomorphism % : 
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-Ko(Vk) — > R- For varieties in Vc, the topological Euler characteristic is an additive 
invariant, and the properties of being additive under decompositions X = YD (X \ 
Y) and multiplicative under products are its defining properties, hence one can 
consider the class [X] in the Grothendieck ring as a universal Euler characteristic 
of the variety X, see [8]. 

Recall that for every prime number p and every n > 1 there is a unique field 
extension F q of degree n of F p , with q = p n . In the case of varieties X defined 
over a finite field F p . One denotes then by X(¥ pm ) the set of F p ™-points of X, for 
m > 1. One then sets 

(2.11) N pm (X) = #X(F pm ), 

the number of F p m-points of X. Then it is not hard to see that N p ™(X) is an 
additive invariant of X, hence it factors through Ko(Vf ), so that we can write 
N p ™([X]) as a function of the class [X] in the Grothendieck ring. 

In the following, we will be considering varieties X defined over Z. In this case, 
we can either regard them as complex varieties by embedding Q C C, or reduce 
modulo p, at the various primes, and obtain varieties X p defined over Fp, for which 
one can consider the counting of points N p ™(X). 

We are interested in the behavior of the number of points N q (X) as a function 
of q. We recall the following terminology (see [7], [35], [35]). 

Definition 2.3. Let X be a variety over Z with reduction X q over F q , for q = p m . 
Then X is polynomially countable if N q (X q ) is a polynomial function of q. 

Because the counting of points over finite fields is an additive invariant, which 
factors through the Grothendieck ring of varieties, one can use polynomial count- 
ability as a test for the motivic nature of the hypersurfaces, as follows. 

If L = [A 1 ] is the Lefschetz motive in K (V), then Z[L] C K (V) is the subring of 
the Tate motives. More precisely, one should view the Tate motives as Z[L, L -1 ] C 
K (M), for M. the category of pure Chow motives, but for our purposes it suffices 
to work with K (V). There is, anyway, a homomorphism Xmot ■ Kq(V) — > K (M) 
defined by the additive invariant of 22 . For varieties that are not smooth and 
projective, one leaves the category of pure motives and one needs to consider objects 
in a more complicated (triangulated) category of mixed motives, which contains a 
(triangulated) subcategory of mixed Tate motives. However, even varieties that 
are singular still define classes in the Grothendieck ring, and the Z[L] part of the 
Grothendieck ring contains the classes of the varieties whose motive is a mixed Tate 
motive. 

Then, if a variety X defined over Z is a mixed Tate motive, it has class [X] € Z[L] 
in Kq(Vz,). Away from up to finitely many primes p (where some bad reduction 
phenomenon can occur), it then follows that N p -m,(X p ) is a polynomial in p m , by 
the fact that the counting N p ™ (X p ) factors through the Grothendieck ring Kq(Vw ) 
and that ^^(A 1 ) = p m . Thus, if a variety X does not have the property of being 
polynomially countable at all but finitely many primes, it follows that it is not a 
mixed Tate motive. 

In fact, polynomial countability at all but finitely many primes is conjecturally 
equivalent to the motive of the variety being mixed Tate: the Tate conjecture 
predicts that determining N p (X p ) for almost all primes p would determine the 
motive of X, see [5J. 
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3. The fibration condition 

In [I] a specific condition for the classes in the Grothendieck ring of the Potts 
model hypersurfaces is identified, called the "fibration condition", according to 
which the the Grothendieck class {Zq} of the complement of the hypersurface Zq 
behaves as one would expect in the case of a fibration on the locus q 0,1. This 
means that the class {Zc, q } is independent of q, for all q ^ 0,1, and that 

(3.1) {Z G } = (T - 1) • {Z GiM } + T# E ^ G K 

This condition is satisfied by the families of graphs considered in [T] , but the gen- 
eral question of whether it holds for more general graphs was not addressed in that 
paper, nor was the question of whether, when the condition on the Grothendieck 
classes is satisfied, the variety Zq is really a locally trivial fibration over the locus 
9^0,1. 

Here we show that, in fact, very simple examples of graphs already do not satisfy 
the fibration condition and the number of points over F p of Zc,q has a nontrivial 
dependence on q ^ 0, 1. 

Proposition 3.1. The graph KA does not satisfy the fibration condition. 

Proof. One can see this by a direct computation of the number of points over F p 
for different values of q. The following table illustrates the failure of the fibration 
condition. 



p 


q 




11 





1771561 


11 


l 


771561 


11 


2 


173799 


11 


3 


173183 


11 


4 


173821 


11 


5 


173513 


11 


6 


174151 


11 


7 


173227 


11 


8 


173447 


11 


9 


173579 


11 


10 


173799 



□ 



By comparison, polygons are the simplest example of graphs that do satisfy the 
fibration condition of pQ. In fact, by Proposition 5.2 of [T] we know that for a 
polygon C m +i with m + 1 sides, one has 

(T — 1 V" — C— 1 V m 

(3.2) {Z Cn+uq ^ tl } = T m+1 + T(T m — (T — l) m ) + ^ ' [ ' , 



with T = L — 1, and it is therefore independent of q =/= 0,1. Comparing this 
expression with Proposition 5.1 of jT], one sees that indeed (3.1) holds for the 
polygons. 
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4. Monte Carlo method for counting points 



Monte Carlo methods for counting points of varieties over finite fields were in- 
troduced in [23] and [57]. We employ here the algorithm described in [33J and we 
show how it compares with the deterministic counting of points, on a sufficiently 
simple example of a graph where both computations can be performed. We use 
again the example of the tetrahedron graph K^, which we already discussed in the 
previous section. 

The following tables provides an explicit comparison between the determinis- 
tic counting and the Monte Carlo method (after 10000, 40000, and 100000 trials, 
respectively) for K4, in the case q = 2. 



p 


Monte Carlo 


% error 


% error bound 


3 


413.9262 


0.002242615 


0.03363242 


5 


4507.8125 


0.013219263 


0.047440384 


7 


20670.9293 


-0.011007641 


0.060059259 


11 


179459.1293 


-0.004845872 


0.079343222 


13 


396763.6998 


-0.015535761 


0.087605918 


17 


1469977.952 


-0.015716002 


0.101770275 


19 


2399339.931 


-0.070218248 


0.108087334 


23 


6928079.605 


0.045288885 


0.119636754 


3 


412.1037 


-0.002170218 


0.01681621 


5 


4453.125 


0.000927175 


0.023720192 


7 


21141.5253 


0.011507837 


0.03002963 


11 


182825.0952 


0.013819407 


0.039671611 


13 


404848.6049 


0.004524793 


0.043802959 


17 


1503167.109 


0.006507159 


0.050885138 


19 


2560472.073 


-0.007777023 


0.054043667 


23 


6661615.005 


0.005085466 


0.059818377 


3 


413.74395 


0.001801332 


0.010635505 


5 


4477.5 


0.006405934 


0.015001967 


7 


20710.92996 


-0.009093825 


0.018992405 


11 


183480.5728 


0.017454225 


0.02509053 


13 


402024.9216 


-0.00248143 


0.027703424 


17 


1499908.538 


0.004325248 


0.032182587 


19 


2532009.315 


-0.018806787 


0.034180216 


23 


6580195.266 


-0.007198912 


0.037832464 



The values should be compared with the deterministic values reported in the 
table in Proposition 5.8 The expected error bound is computed as in [23], with 
an error margin e of y/4b(log(2/S))/N, where b is the fraction of points that are 
roots, N is the number of trials, and 1 — <5 is the probability of being within e of 
the correct fraction. 

The Monte Carlo method is useful to test properties such as the fibration con- 
dition, when direct deterministic computations become intractable. In fact, when 
using the Monte Carlo method in cases with non-constant g-dependence, the de- 
viance from a constant value is typically greater than what allowed by the error 
estimate. 
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5. Failures of polynomial countability 

We now analyze failures of polynomial countability for the Potts model hyper- 
surfaces, which confirms the fact that these hypersurfaces are much less likely to be 
mixed Tate motives than the graph hypersurfaces of quantum field theory, although 
families of mixed Tate cases can still be constructed, using the deletion-contraction 
relation derived in [1], in cases where a good recursive procedure is possible, such 
as the polymer chains analyzed in pQ. 

We use a variant of the algorithm used by Stembridge in the case of the 
graph hypersurfaces. This can be adapted to our setting, because some general 
results proved in [36] apply both to the graph polynomials and the multivariate 
Tutte polynomials, in particular those discussed in Lemma |5.6| below. 

5.1. The special case q = 1. There is a very special case, which can be treated 
separately by direct geometric reasoning, which is the case with q = 1. Notice that 
this is not a physically relevant case, since q is the number of spin states. 

When q = 1, the multivariate Tutte polynomial assumes the very simple form 

z G ( q ,t)\ q=1 = n*«= n( i+t ^- 

A(ZE e£A e£E 

In fact, multiplying this out gives all possible combinations of the i e 's, which cor- 
respond to all possible subsets A of E. 

Lemma 5.1. In the case q = 1, the counting of points over F p is given by 

N p (Z G . q=1 )^p# E ^-(p-l)# E ^. 

Proof. The edge weight tuple satisfies the Tutte polynomial iff at least one of the 
edge weights t e is —1. Therefore, we can get the number of roots by taking the 
total number of tuples and subtracting the number of tuples where none of the edge 
weights is —1. If there are p elements in the field and m edges, the number of roots 
is p m — (p — l) m . Therefore, the number of roots always depends polynomially on 
p when q = 1. □ 



One can also see this directly from the class [Zq q=1 ] in the Grothendieck ring, 
as observed in pQ, the locus of zeros of the polynomial Ilee£;(l + is isomorphic 
to the union of coordinate hyperplanes, whose complement in A# B ( G ' is a torus, 
whose class is T# E ( G \ Thus, the class of the complement is {Zc. q =i} = T# E (°) 
and [Zq^i] = \J^ E ^ — T# E ( G \ which immediately gives back the formula of 



Lemma 5.1 since N v (h) = p and -/V p (T) 



5.2. A trivial failure for q ^ 1. We now consider the case where q ^ 1 (we also 
assume q ^ 0). In this case, we show that there is always a "trivial" failure of 
polynomial countability, which has to do with the primes that divide q. In fact, in 
the way it is originally written, the multivariate Tutte polynomial has an overall 
multiplicative factor of q kl - G \ with k{G) the number of connected components of 
G. Thus, whenever p is a prime that divides q, the polynomial becomes identically 
zero in the reduction mod p. 

We start with the following simple observation. 
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Lemma 5.2. // a set of points satisfies a single variable polynomial at infinitely 
many values but does not satisfy that polynomial at certain other values, the entire 
set of points cannot satisfy a single variable polynomial ( even if it is a different 
polynomial) . 

Proof. Suppose given an infinite set S with an infinite subset T C S, such that 
the points of T satisfy a polynomial Q(x) of degree m and the points of S satisfy 
a polynomial P of degree n, then these two polynomials must be the same. In 
fact, both polynomials have to agree on the infinite set T, hence their difference 
R(x) — P(x) — Q(x) is a polynomial with infinitely many zeros, that is, the trivial 
polynomial R = 0. □ 

We then have the following general behavior. 

Corollary 5.3. For a fixed value q =/= 0, 1 of the number of spin states, let G be a 
graph for which the class [Za. q ] in Kq(Vz) is in Z[L]. Then the counting function 
N p {Z G>q ) is given by 



(5.1) N p {Z Gs ) = 



p#£(G) p\q 

P q (p) P k 



where \Z G ^ q \ = P g (L) — oq + a{L + • • • + a^h N in Z[L]. Thus, N v (Za, q ) is not a 
polynomial function. 

Proof. If q = in the field (that is, if p is a divisor of q), everything is a root, so 
there are p m roots where m is the number of edges. If q (p does not divide q), by 
the form \Z Gq ] = P q (h) = oq + a{L + • • ■ + a^rL^ of the class in the Grothendieck 
ring we obtain that Np{Za, q ) — Pq(p)- The last observation then follows from 
Lemma 15.21 □ 



Notice that this type of failure of the polynomial countability, due only to the 
primes that divide the number q of spin states, should be regarded as an "acciden- 
tal" and not a "serious" failure of polynomial countability, in the sense that, for 
example, it does not really affect the nature of the classes in the Grothendieck ring. 

We can see this explicitly in some simple examples. Consider first the case where 
G = T is a tree. Then the multivariate Tutte polynomial is of the form 

(5.2) Z T (q,t)=q JJ (q + t e ). 

e£E(T) 

Example 5.4. The counting function Np(Zr miq ) for a tree T m with m edges is 
given by 

(5-3) N p{ Z Tm , q ) = [ j£_ (p _ ir % 

Another very simple explicit example is that of polygons, for which the expression 



(3.2) for the class in the Grothendieck ring gives the following. 



Example 5.5. The counting function Np(Zc m+1 , q ) for a polygon C m +i is given by 

pm+i ^gjj p|g^ while when p J(q, it is given by 

(5.4) 

(b — 9\ m — (— 1 ) m 

N p (Z Cm+uq ) = p"+i-((p-i)"+i + (p-i)((p-i r -(p-2n+ " 1 ] ). 
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Another way to see that N p (Z Gm+lyq ) is a polynomial in p, for all p that do 



not divide q, which does not use directly the expression (3.2) for the class in 
the Grothendieck ring, is to show, as in Stembridge [3B] that a partial deletion- 
contraction relation holds at the level of the counting functions. This reflects the 
partial deletion-contraction relation at the level of classes in the Grothendieck ring 
proved in pQ. 

5.3. Probabilistic counting functions. We recall the relevant notation from 
[35] . Given a set of polynomials {/i , . . . , fk} in Z[x\ , . . . , x m ] , one writes 3 [/i fk] 
for the probability that all the /j vanish at a uniformly randomly chosen (xi, . . . , x m ) G 
Fq. Then, as a function of the field cardinality q = p r , one has 

(5.5) 3[/i,...,,/ fe ](q) = q- m iV q (X p ), 

where X p is the reduction mod p of the variety denned over Z by the polynomials 
{fx, ...,fk} and iVq is the number of its F q -points, for q a power of p. Then 
by Proposition 2.1 of [33], the variety X is polynomially countable if and only 
if 3[/i) • • • , j /fc](q) is a polynomial in q _1 . We will refer to 3[/i, • • • , , fk] as the 
"probabilistic counting function" . 

Then we have the following analog of Lemma 3.2 of [36] . 

Lemma 5.6. Let Z Gtq = Z G (q,t) denote the multivariate Tutte polynomial, for 
fixed q ^ 0, 1, as a function of the edge weights t = (t e ), and let ~5[Z G ^ q ] be the 
corresponding probabilistic counting function. Then this satisfies the relation 

(5.6) 3[Z G , g }(p) = - p-^G/e,,] + 3[^G/e, g! ^Gse,,]- 

Proof. The multivariate Tutte polynomials satisfy a deletion-contraction relation 

(5.7) Z G (q,t) = Z G U<l,i (e) )+teZ G/e (q,i^), 

for any edge e (regardless of whether it is a regular edge, a bridge, or a looping 
edge), with t^ denoting the set of edge weights with t e removed, and where G/e is 
the graph obtained by contracting e in G and G \ e the graph obtained by deleting 
e in G. 

Then the argument of Proposition 2.3 of [36] and the deletion-contraction rela- 



tion (5.7 1 show that one obtains the relation (5.6). □ 



A more refined deletion-contraction relation of a similar nature to (5.6) was 
proved by algebro-geometric methods in [I], for classes in the Grothendieck ring. 
Namely, for the classes of the hypersurface complements one has 

(5-8) {Z G<q } = L • {Z G/e , q n Z G ^ e , q ] - {Z G /e,g}- 



One sees then that one can also recover ( 5.6 ) directly from ( 5.8 1 . Indeed, we have 



from (5.8), that (away from primes that divide q) the counting function satisfies 

N p (A# E ( G ) x Z G . q )) = pN p (A* E ( G ^ x {Z G/e , q n Z G ^ q ))) 
- ^(A# £ («)-K2 G/ei! )), 

which gives 

N p (Z G , q )) = pN p (Z G/e , q n Z G ^ q ) + p# £ ( G )" 1 - N p (Z G/e , q ). 
Thus, one obtains, as expected, 

3[ZgJ(p) = p-# E(G) N p (z G , q ) = p- 1 + 3[z G/e , q n z G ^ q ] - p-^iZc/eJ- 
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Similarly, one can define another probabilistic counting function, which we de- 
note here by 3 v [/i, • • • ,/fc] = 1 — 3[/i, ■ • ■ , fk]- This is the probability of not 
having a common zero of the /j. Just as 3[/i, ■ • ■ i/fcKl) i s related to the num- 
ber of points JVq(Xp) by 3[/i, • • ■ , /fc](q) = c\~ m N q (X p ), the probabilistic counting 
function 3 v [/i, • • ■ , fk] satisfies 

(5.9) 3 v [/i, ■ ■ • , /*](q) - q" ro (q m - N q (X p )) = q~ m N q (A m v X p ). 

In other words, 3 v [/i, ■ ■ • , fk] is the probabilistic counting function of the variety's 
complement. Thus, from (5.8) we have 



(5.10) y[Z G , q ] = y[Z G ^Z G/e , q ] - q- 1 3 v [^ G / 



5.4. The case of the graph K4. We now focus, in particular, on the polynomial 
countability question for the case of the tetrahedron graph K4, and we show that 
it exhibits a more serious failure of polynomial countability, which is not only due 
to the primes that divide q. 

Lemma 5.7. "$\x 2 + 2x + 2] is not a polynomial. 

Proof. First observe that there are infinitely many primes of the form p = Ak + 
3. This follows from Dirichlet's theorem, showing more generally that there are 
infinitely many primes in any arithmetic progression. (For this particular case 
there is also a direct elementary proof, which we do not report here.) Then notice 
that, for all primes p of the form p = 4fc + 3, the polynomial x 2 + 2x + 2 has no 
solutions over F p . In fact, by the quadratic formula, x 2 + 2x + 2 has a solution if and 
only if x = ~ 2± v / 4-8 _ — 1±-\/ — 1. This has a solution iff G F p . Assume this is 
true, and that p = 4fc + 3. Let a = then we have a 13 " 1 = a 4k+2 = (a 2 ) 2k+1 = 

_ 1 2fc+i = _ 1) but by p ermat ' s Li tt i e Theorem, a^ 1 = 1. Therefore, x 2 + 2x + 2 
has no solution for all primes of the form p = 4k + 3. Now, if 3[# 2 + 2x + 2] were 
a polynomial, then it would have infinitely many zeros, since there are infinitely 
many primes p (those of the form p = 4fc + 3), where x 2 + 2x + 2 has no solution, 
but })[x 2 + 2x + 2] is not identically zero: for instance, when p = 5, there are two 
roots, x — 1 and x = 2. Thus, ^>[x 2 + 2x + 2] is not polynomial. □ 

Consider then the graph KA and in the case with q = 2 (Ising model). A direct 
calculation based on a version of the algorithm of Stembridge gives 
(5.11) 

N p (Z Kit9= a) = p 5 3[P(z)] + 2p 3 3[2 + 2x 4 + x 2 ] 

+2p 3 3[2 + 2x 5 + a:§] + p 5 - p 4 - 3p 3 + 13p 2 - p - 1 
-p 5 3[2, P{x)] - 2p 3 3[2, 2 + 2x 4 + x 2 ] - 2p 3 3[2, 2 + 2x 5 + x 2 ] 
+3[2](p 6 - p 5 + p 4 + 3p 3 - 13p 2 + p + 1), 
where the polynomial P{x) is of the form 

P(x) — 4X2X4X5 + 8X2X3X4X5 + 4x2X3X4 + 4X2X3X5 + 4x3X4X5 

+4x3X4 + 4x2X5 + 2x 3 x 4 + 2 

+X2X3X4 + 2x2x3X4 + 2x3X4X5 + 2X3X4X5 + X3X4X5 
+2X3X4X5X3 + 2x2X4X5x3 + X3X4X5X2 + 2X3X4X5X2 

+2x3X5 + 2x3X4X5 + 2x3X5X3 + 2x 2 x 2 + XjXg + X2X4X 2 
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X2X4X5X3 + 2X2X4X^X3. 



+2X2X4X 5 + X 2 XsX 5 + 2X2X3X 5 

The 3 [2] terms are only 1 in the case where 2 = 0, which we can ignore as we are 
looking for non-trivial failures of polynomial countability. Among the remaining 
terms, we have seen in Lemma 5.7 that the 3[2 + 2x + x 2 } contribute non-polynomial 
expressions. This strongly suggests that 3[^.K"4,g] itself may be non-polynomial. 
However, one needs to make sure that there are no cancellations (however unlikely) 
coming from the term p 5 3[-P (x)). 

This is confirmed also by looking directly at the number of points for sufficiently 
many primes, as follows. 

Proposition 5.8. For q = 2 (Ising model case), the function Np(Zx4. q =2) is non- 
polynomial. 

Proof. The values for the first few primes p ^ 2 give the following table. 



p 


3 


5 


7 


11 


13 


17 


19 


23 




413 


4449 


20901 


180333 


403025 


1493449 


2580541 


6627909 



p 2 a 2 + p 3 a 3 + p 4 a 4 + p 5 a 5 = N p , 



Proposition 2.2 of [36 shows that, since the graph K4 has six edges, if N p (Z K4: , 9=2 ) 
is given by a polynomial, then this polynomial would have degree at most five. Thus, 
one obtains from the table above a system of linear equations for the coefficients of 
this polynomial, 

(5.12) a + pdi 

for p and Np in the table above. Already when solving the first five equations in 
this system, one finds only a solution with rational, non-integer values of the dj. By 
Proposition 6.1 of [31] . one knows that if Np(ZK4, q =2) is polynomial with rational 
coefficients then the coefficients must in fact be integers. Indeed, if one further 
considers the remaining three equations from the table of values above, one finds 
that there are no solutions. Thus, Np(ZK4, q =2) is non-polynomial. □ 

While the result of Proposition |5.8| alone suffices to show that the Potts model 
hypersurface of the graph KA is not polynomially countable, we have included 
the previous discussion to illustrate a possible explicitly source of non-polynomial 
terms. 

5.5. Normalized Tutte polynomial. We make a small additional remark on the 
trivial failures of polynomial countability. In order to avoid the presence of an 
overall factor g fc ( G ) in the multivariate Tutte polynomial, one sometimes considers, 



instead of the polynomial (2.4 1, the normalized version 

(5.13) Z G (q,t) = q- k ^Z G (q,t). 

The effect of this change on the counting function is only to alter it at the primes 
p that divide the number of states q. For example, in the two simple examples of 
C4 and K4, one finds the following. 



Example 5.9. For the normalized polynomial (5.13), the counting function for C 4 
with q — 2 is given by 



P 



3p 2 3[2]+5p 2 + 5p3[2]-7p-l, 



which differs from the one for the Tutte polynomial (2.4) by 

3[2](~p 4 +p 3 + 2p 2 -2p-l). 
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For the graph K4 with q = 2 the counting function for the normalized case gives 
p 5 3[P{x)} + 2p 3 3[2 + 2x4 + x\] + 2p 3 3[xl + 2x 5 + 2] 
+p 5 _ p 4 _ 3p 3 + Up 2 _ p _ 1 + 3[2](2p 3 - 15p 2 + 2p), 

which again differs from the case of the original Tutte polynomial only at p = 2, 
with the difference given by 

3[2](-p 6 +p 5 -p 4 -p 3 - 2p 2 +p- 1) +p 5 3[2,P(x)] 

+2p 3 3[2, x\ + 2x 4 + 2] + 2p 3 3[2, x\ + 2x 5 + 2}. 

6. Thermodynamic averages and periods 

Coming back to the physical interpretation of the multivariate Tutte polynomial 
Zg(q> t) as the partition function of the Potts model with q spin states on the graph 
G, the edge weights are of the form t e = exp(/3J e ) — 1, where J e > is the energy 
(in the ferromagnetic case) and (3 is the inverse temperature. One knows that, for 
a finite graph and in the ferromagnetic case, there are no phase transitions, that 
is, no zeros of Zc{q,t) that fall in the range t e > and the phase transitions only 
appear in the limit n — > oo of a family of finite graphs G n approximating an infinite 
graph (see [53] )• 

Given a function f(J) = F(t) of the energies J = (J e ), or of the edge weights 
t = (t e ), one can compute the thermodynamic average 

((\\\ !F\ - £AC£g fe(A) * , fcOII e gA*e _ 1 V k(A) p/, \ TT ■ 

((U) <F> - w^u; {A) lk 

where we write F(tA) = F(t)\ ti!=0 y e ^ A . We are especially interested here in the case 
where F(t) is a polynomial function of the edge variables, with rational coefficients. 

Moreover, we can further average over a range of energies (at a fixed tempera- 
ture), for example, by letting the edge variables range over the simplex A = {t = 
(te) | te > 0, 2^ e te = !}• One then finds an expression that is formally very similar 
to the parametric Feynman integral in perturbative quantum field theory, namely 

(6-2) TXT / (F) dv = [ dv(t), 

V ; Vol (A) J & 1 Vol(A)J A Z G (q,t) u ' 

where we write Pg.fC?^) f° r the polynomial 

ACE eeA 

under the assumption that F is itself a polynomial. 



The normalization factor in ( 6.1 ) can be computed easily: the volume of a regular 



n-dimensional simplex with side length a is given by the well known expression 



Val(A n (a)) = ^V^- 

Thus, the interesting number, whose nature one wants to investigate, is the remain- 
ing expression 

r - 

A Z G (q,t) 



(6.3) / 
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Figure 1. A polygon polymer chain (from pQ). 



We see that the expression (6.3) that we obtain in this way is formally similar 



to the parametric Feynman integral (although with a very different physical inter- 
pretations), with Zc{q 1 t) the multivariate Tutte polynomial instead of the graph 
polynomial, and with Pg,f instead of the second Symanzik polynomial. 

An interesting difference with respect to the quantum field theory case, which 
makes the Potts model case nicer, is that (as recalled above) the partition function 
Zc(q,t), for a finite graph G does not have any zeros in the domain of integration. 
Thus, these integrals are genuine periods and are not plagued by the infrared di- 
vergences problem that occurs in their quantum field theoretic analogs, where the 
graph hypersurface intersects the simplex A. 

By analogy with the context of quantum field theory, one can then ask the 
question of what kind of numbers one obtains by evaluating the integrals (6.3|, for 
arbitrary polynomials F(t) with rational coefficients. 

The integral (6.3) is a period of the variety A™ \ Zc, q , where n — #E(G): the 



integral of an algebraic differential form over a domain of integration given by a 
(semi) algebraic set defined by algebraic equations. Periods are a very interesting 
class of numbers, see the detailed account given in [28] , 

In particular, what kind of numbers can arise as periods of a given variety de- 
pends on the nature of the motive of the variety. For example, it was long conjec- 
tured and recently proved in |13j that the periods of mixed Tate motives over Z 



of the form 



-linear combinations of multiple zeta values. The latter are expressions 



-) 



E 

0<fci<...<fe r 



> 2. 



with integers ni > 1 and n r 

Thus, for example, in the cases of polygon polymer chains considered in [T], 
where it is shown by the explicit formula (3.2 ) that the motives involved are mixed 
Tate, one can conclude that the periods (6.3) must evaluate to combinations of 



multiple zeta values, which is not obvious from the definition of the integral. 

We assume here some standard conjectures about motives, in particular the Tate 
conjecture which implies that determining N p (X) for almost all primes p determines 
the motive of X (see [B], §1.1). 



Proposition 6.1. Let ^^Gjq be a chain of linked polygons, obtained by joining N 
polygons, each with m + 1 sides and connected by a chain of k egdes, as illustrated in 
Figure [7J Then the thermodynamic averages (6.3) are 
of multiple zeta values. 



-linear combinations 
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Proof. The explicit formula given in Proposition 5.5 of pQ for the class in the 
Grothendieck ring gives 



where {Zc m+lt q} is as in (3.2 1. This shows that the class in the Grothendieck ring 
of the variety Z{ m ,h)Q N „ is in the mixed Tate subring Z[L] C Kq(Vj,). In principle, 
this does not yet say that the motive of Z( m ,k)Q N „ is mixed Tate (in the category 
of mixed motives), as the class [Z( m ,k) GN q ] determines only its motivic Euler char- 
acteristic Xmot(Z(.m,k)Q Ntq ). However, if the conjecture holds, which predicts that 
knowledge of the counting function Nn suffices to determine the motive, then we 
can conclude from the expression above for the Grothendieck class that the motive 



is mixed Tate. Then the main result of [12] shows that the period (6.3) has to be 



a combination of multiple zeta values as stated. □ 

One can then ask a more precise question. 
Question 6.2. What combinations of multiple zeta values arise in the period com- 



putations (6.3) for the polymer chain graphs ^ m ' fc ^Gjv? How do they behave as the 



parameters m,k,N get large? 

With respect to this last question, and the more general question of identifying 
criteria for when a graph satisfies the polynomial countability condition, it is worth 
mentioning another important difference between the quantum field theory case 
and the Potts model hypersurfaces. 

In the case of the graph hypersurfaces of quantum field theory, where the hy- 
persurface is defined by the vanishing of the Kirchhoff polynomial, one can use the 
matrix-tree theorem to obtain some more precise information on the periods. In 
fact, the Kirchhoff polynomial can be written as ^g = det(Mct), where A4q IS 
a matrix defined in terms of the combinatorics (the incidence data) of the graph. 
One then has algebraic relations between determinants of minors of the matrix Mg, 
which define relations between the corresponding Dodgson polynomials ^q J of the 
graph. These identities were used in |14) to identify structural properties of graphs 
whose associated Feynman period is a multiple zeta value. This identifies a "linear 
rcducibility criterion" for graphs, which implies the mixed- Tate property. Another 
setting in which the matrix-tree theorem for Kirchhoff polynomial can be useful in 
the quantum field theory case is the result of [5J, where one uses the determinant 
description of the Kirchhoff polynomial to map the Feynman integral computation 
to the complement of a determinant hypersurface. In the case of the multivariate 
Tutte polynomials, however, it is well known that one does not have an analog of 
the matrix-tree theorem. In fact, the problem of computing the Tutte polynomials 
is known to be #P-hard, see [55] . 

6.1. Potts models on tetrahedral chains. Our previous observations on the 
tetrahedron graph suggest that the same approach we just described for polygonal 
polymer chains may lead to possibly very different answers when applied to a dif- 
ferent type of chain graphs, the tetrahedral chains that arise naturally in the study 
of silicates. 

The simplest example of Potts model, the Ising model with q = 2 spin states, has 
been extensively used in studying the properties of tetrahedral chains that arise in 
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FIGURE 2. Tetrahedra in a single-chain configuration. 



silicates and Si-Al minerals (see §6 of [3U]). Example of such Ising model analysis 
can be found in [TJ]], [2"T] . 

Here we consider the type of single chain tetrahedral configuration, as illustrated 
in Figure [2j This is realized, for instance, by inosilicates such as (Si03~) n . 

Recall that, for a graph G that is obtained as a chain of graphs G\, . . . , Gn, where 
Gi and G^+i have a single vertex in common, the classes in the Grothendieck ring 
of the corresponding Potts model hypersurface complements (with q ^ 0) is simply 
a product (see Lemma 3.5 of pQ) 

N 

(6-4) {Z G J = Y[{Z Gi , q }. 

i=i 

Thus, if Gjv is the graph obtained by linking together N copies of the tetrahedron 
graph K4, as in Figure [2j the corresponding class for q ^ is given by {Zc N , q } — 
{Zk4, q} N , or equivalently 

[Z GN , q ]=h m -(h 6 -[Z KA J) N , 

which implies that the counting function satisfies 

N p (Z GN , q ) = p 6N - (p 6 - N p (Z K4 , q )) N , 

at primes p that do not divide q. In the case of the Ising model, where q = 2, the 



counting function Np(ZK4,,q=2) is the one computed in (5.11). 

Thus, one sees that, unlike the case of the polymer chains analyzed in 1], that 
involved the polynomial counting functions of polygons, here the counting function 



involves non-polynomial expressions, as shown in Proposition 5.8 Although this 



does not prove that the motive is non-mixed Tate, as there could be cancellations 



coming from the other terms in (5.11), it appears to be a good candidate to test. 



So one can ask the following question. 

Question 6.3. Let Gn be the graph obtained by linking together N tetrahedra as 
in Figure [M Are there choices of polynomials F(t) with rational coefficients, for 



which the period (6.3), computed for this graph is not a Q[-^ 7 ]-linear combinations 



of multiple zeta values ? What kind of period is it, if that is the case ? 

Notice that the analogous question in the quantum field theory setting is not yet 
solved, even though there are explicit non-mixed Tate examples of graph hypersur- 
faces as in [16] , In fact, there is not yet an explicit example of graph hypersurface for 
which the Feynman amplitude itself would be a non-mixed Tate period, although 
there are good reasons to expect that, as explained in [Tt)] . 
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